Application of the Fròbenius method 

to the Schròdinger equation for a spherically 

symmetric potential: anharmonic oscillator 

Przemystaw Koscik and Anna Okopinska 
Institute of Physics, Pedagogical University, 
Swietokrzyska 15, 25-406 Kielce, Poland 

E-mail: koscik@pu.kielce.pl 
Abstract. 

The power series method has been adapted to compute the spectrum of the Schròdinger 
equation for central potential of the form V(r) = -0- H — — + dir 1 . The bound- 
state energies are given as zeros of a calculable function, if the potential is confined 
in a sphcrical box. For an unconfincd potential the interval bounding the energy 
eigenvalues can be determined in a similar way with an arbitrarily chosen precision. 
The very accurate results for various spherically symmetric anharmonic potentials are 
presented. 
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1. Introduction 

The exact solution of the Schròdinger equation can be obtained only for a few particular 
forms of potentials, in other cases one has to resort to approximations or numerical 
techniques. Many approximation methods have been developed for solving problems 
in one-dimensional space. Approximate solutions to the Schròdinger equation have 
been also studied for spherically symmetric potentials in D— dimensional space, both 
by methods elaborated for one-dimensional space, e.g. the Hill determinant method 
[T], the variational approach |2j, and by methods dedicated to D-dimensional problems, 
e.g. the shifted 1/D expansion [SUI]. Here we show that highly accurate solutions to 
the Schròdinger equation can be determined for various types of spherically symmetric 
potentials with the use of the Fróbenius method (FM). The method consists in expanding 
the solution of a differential equation into power series 0, and was originally applied by 
Barakat and Rosner jü] to compute the spectrum of one-dimensional quartic oscillator 
confined by impenetrable walls at x — ±R. The energy eigenvalues of the system 
have been obtained numerically as zeros of a function, calculated from its power series 
representation. Moreover, it has been shown that the bound-state energies of the 
confined system approach rapidly those of the unconfined oscillator for increasing R. 
Low-lying eigenvalues for other one-dimensional potentials j?j have been also successfully 
calculated in a similar way Recently, a modified treatment of unconfined systems 
allowed for a very accurate determination of the ground-state energy for the quartic 
oscillator [B] . In all the cases studied the potential was a finite function, and a solution 
was expanded around an ordinary point of the differential equation. Here we study the 
application of the FM for solving the radial Schròdinger equation, which requires that 
an expansion around a regular singular point should be used. 

The outline of the present work is as follows. In Section|2]the solution of the radial 
Schròdinger equation in form of a generalized power series is discussed. The case of a 
spherically symmetric potential bounded by an impenetrable wall at r = R is studied in 
SectionEl In this case, the energy eigenvalues can be easily determined by finding the 
roots of the polynomial, which is illustrated on the example of the confined harmònic 
and anharmonic oscillators and Hulthén potential. The case of unconfined system is 
studied in SectionEJ where a scheme for determining an arbitrarily large set of bound- 
state energies is developed. After demonstrating the performance of the method on 
the exactly solvable example of the Kratzer potential, the results for the unconfined 
oscillator are presented for various choices of anharmonic parameters. 

2. Expansion around a regular-singular point 

The Schròdinger equation for a spherically symmetric potential in 3— dimensional space 
can be reduced to an ordinary differential equation in the radial variable 




(1) 
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where l is the angular momentum quantum number, and the units h = 1, m = 1 are 
used. Upon introducing the function u(r) = rR(r), the differential equation (0) takes 
the form of the one-dimensional Schròdinger eigenvalue problem 



-n—2+Veff(r,l) 



u(r) = Xu(r), (2) 



where the effective potential reads 

V eff (r,l) = l -^ + V(r). (3) 

The point r = is a regular singular point of the radial equation, if the potential 
V(r) diverges but r 2 V(r) remains finite as r — > 0, which is the case for the interaction 
potential of the form 

V(r) = ^ + d -f + V re9 (r), (4) 
where the regular part is represented by a convergent series 

oo 

U0 = E^ ( 5 ) 

i=0 

In this case the FM can be applied with the radial wave function represented as a 
generalized power series 

oo 

u(r) = r 5 E a i r \ (6) 

i=0 

where a ^ 0. In what follows, we will take üq = 1, since normalization of wave function 
is irrelevant in our calculation. Substituting (jüj) into (J2J) we obtain the equation 

-1 oo oo oo oo 

-rE[( i + í )( í + 5 - 1 )-Kí+l)]«^ í +(E^)(E a ^) = AE^ í+2 >(7) 

^ i=0 i=0 i=0 i=0 

which, by comparing the coefficients of like powers of r, yields the recurrence relation 

i 

[{i + 8){i + 8-l)-l{l+ 1)]oí - 2 E ^-2- n a„ + 2a,_ 2 A = 0, (8) 

n=0 

where <2j = for i < 0. Setting z = in the above relation, we obtain the indicial 
equation 

[5(6 - 1) - 1(1 + 1) - 2d_ 2 ]a = 0, (9) 
which is solved by 

6, = l -(l - ^8d_ 2 + (l + 2/) 2 ), 5 2 = i(l + v /8rf_ 2 + (l + 2/) 2 ). (10) 

The Fuchs's theorem jHj asserts that the generalized series © converges, and 
both linearly independent solutions of the Schròdinger equation (J2J are obtained as 
generalized series, one with 8 — 8\ and the other with 5 = 62 (except the special case 



when 82 — 81 = ySd_2 + (1 + 2Z) 2 is equal to a non-negative integer). The value of 8 
determines the behavior of u(r) for r — > 0, and only 8 > ~ is acceptable j^j, since only in 
this case the mean value of the kinetic energy is finite. Such a solution to the Schròdinger 
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equation (J2J) exists only if the potential is such that d_ 2 > — §• This solution contains 
only the series with 6 = 62, and in the following will be denoted by 



30 



u(r, A) = / 2 2 ai r\ (11) 

where the dependence on the energy eigenvalue À is explicitly marked. The coefncients 
of the series are determined recurrently as 

a _ ~ 2ÀOi_2 + 2 5Zn=0 di-2-nO-n /.^N 

i(i + ij8d_ 2 + (l + 2l) 2 ) 

The function u(r, X) has an important property that u(0, A) = 0, which can be regarded 
as a boundary condition for the radial equation at r = 0. The second boundary condition 
should be chosen according to the physical bounds in the investigated problem. 



3. Confined potentials 

The simplest scheme arises in case of radially symmetric potential V(r) additionally 
bounded by an infmitely high wall at r = R. In this case the second boundary condition 
for a particle with angular momentum / is of the form 

u(R,\) = 1 (13) 

yielding an exact quantization condition that gives the bound state energies as zeros of 
a calculable function. The vàlues of À, determined from the above condition, are will be 
denoted by X n i, where n = 0, 1, ... counts the number of zeros in the radial variable in 
the Sturm-Liouville eigenvalue problem jTÜ]. In the numerical calculation, the function 
u(R, X) has to be approximated by truncating the series in (fTT]) at suitably high order K, 
which can be done with an arbitrary accuracy, as the series is convergent. The truncated 
function is a polynomial of degree K in the variable A, that can be used to advantage 
in determining the numerical vàlues of its zeros. In this way an increasing amount 
of bound-state energies can be obtained with increasing accuracy as the truncation 
order K —>■ 00. Considering the growing interest in quantum confined systems, the 
method (and its possible modifications for different symmetries of the confining box) 
can find application in the studies of semiconductor nanostructures such as quantum 
dots El E] • Here we consider a few examples. 



3.1. Harmònic Oscillator 

The performance of the FM will be demonstrated for the spherical harmònic oscillator 

V(r) = yr 2 , (14) 

enclosed by a sphere of the radius R. The typical behavior of u(R, X) as a function 
of À is shown in FigC] in the example of the oscillator with tv — 1 for l — and 
R = 2.5. The approximate u(R, X) is a polynomial in A and the truncation affects only 
its behavior for large A. We determine thus the bound state energies as roots of the 
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obtained polynomial with the use of NSolve procedure from the Mathematica package, 
without any need for introducing starting vàlues for A. The stability of the numerical 
results was achieved by increasing the number of non-vanishing terms K until the vàlues 
of X n i(R), corresponding to the states (n,l) of angular momentum l, become stable to 
the desired accuracy. The vàlues of X n i(R), obtained for the harmònic oscillator of 
frequency ou — 1, are compared in Table ^ with the eigenvalues of the unconfmed 
oscillator, E ni — 2n + l + |. In FigO the low-lying states energies are plotted as a 
function of the confmement radius R . One can observe how the non-degenerate levels 
of the confined system approach the equidistant states of the unconfmed oscillator, as 
the radius of enclosure grows. 




uH2 . 5, 1L 



2 4 6 8 10 12 

1 



Figure 1. Behavior of u(R, A) for the spherically symmetric harmònic oscillator with 
w = 1 and l = at R = 2.5. 



Table 1. The low-lying energy levels of the harmònic oscillator with uj = \ confined 
at various R. 



1 


n 


A n *(l-5) 


A n *(2-5) 


A n z(3) 


•^nz(3.5) 


A ni (4) 


E n l 








2.5049762 


1.5514217 


1.5060815 


1.5003995 


1.5000146 


1.5 


1 





4.9035904 


2.6881440 


2.5312925 


2.5029102 


2.5001438 


2.5 





1 


9.1354221 


4.1842613 


3.6642196 


3.5233023 


3.5016915 


3.5 


2 





7.8717305 


3.9535289 


3.5982477 


3.5125803 


3.5008421 


3.5 



3.2. Anharmonic Oscillator 

The quantum anharmonic oscillator is widely used to describe the physical phenomena, 
especially in the condensed matter and molecular physics. However, the case of 
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Figure 2. Behavior of X n i(R) — Xqq(R) for the harmònic oscillator in function of the 
confinement r adius R. 



a confined anharmonic system has been discussed only for the one-dimensional 
example [§J. Here we study the spherically symmetric anharmonic oscillator with a 
potential of the form 



V AO (r) 



LO 



-r + gr 



2.1 



(15) 



for different vàlues of the power J . The analysis can be simplified by using the rescaled 
variables 



1 ^ 2 

r = g 2J + 2 r, and A = Xg 2J + 2 



then Eq. (J2J) takes a form 



1 d 2 1(1 + 1) ^ 



u r 



Àií(f), 



where 



zf 2 + f 2J 



(16) 



(17) 



depends on the dimensionless parameter 

Z 



LO 4 
2J+2. 



which accounts for a relative strength of the harmonicity and anharmonicity. In the 
following, we skip the hats over r and À. We consider the solution of (|16jh which is 
regular at r = 0, as given by the generalized power series 



u (r, A) = r l+1 air 2i 

with the recurrence relation of the form 

_ -2Aaj_i + 2zai-2 + 2ai_^j +1) 
üi ~ 2z(l + 2i + 2l) 



(19) 



(20) 
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where a = 1. In the case of the oscillator enclosed by a sphere of the radius R the 
bound-state energies are easily determined as zeros of u(R, A). In TableElthe numerical 
results are presented on the example of the quartic oscillator (J = 2) with different 
vàlues of the parameter z and for different radii of enclosure R. The spectrum of an 
unconfmed anharmonic oscillator will be studied in the next section. 



Table 2. The bound-state energies X n i(R) of the confined anharmonic oscillator 



l 


n 


z 




A„,(0.5) 


A„z(l) 


A ní (1.5) 


A ní (2) 


Anz(2.5) 








5 


20 


.09875027 


6.39291817 


5.10864245 


5.06974045 


5.06966368 






4 


20 


.02840681 


6.13228781 


4.69129108 


4.63348473 


4.63330831 






3 


19 


.95799389 


5.86753935 


4.24715028 


4.16026790 


4.15984713 






-1 


19 


.67564622 


4.76592159 


2.12668280 


1.63322595 


1.61278177 






-3 


19 


.53405395 


4.18860246 


0.80447958 


-0.38394286 


-0.54059117 


1 





5 


40. 


.85989845 


12.08775101 


8.80212012 


8.63744701 


8.63688371 






4 


40. 


.76648791 


11.73035216 


8.16458083 


7.93950715 


7.93833601 






3 


40. 


.67302607 


11.36972235 


7.49829458 


7.18963353 


7.18714592 






-1 


40. 


.29866554 


9.89468226 


4.51245222 


3.40693967 


3.34739553 






-3 


40 


.11117731 


9.13754642 


2.81077229 


0.78285838 


0.50413322 



3.3. Hulthén potential 

The FM can be also useful for computing energy eigenvalues, in case when the series 
expansion of its regular part ((SJ) converges only in a finite interval py. In this case, the 
convergence of the generalized series solution ()11|) is granted only for r < py and the 
problem is well defined only for the confinement radius R < py. If the confining box 
projects beyond the convergence sphere, we must be very careful, since for r > py the 
series representat ion u(r, A) does not necessarily coincide with the solution of the radial 
Schròdinger equation in the potential V(r). As an example, we consider the Hulthén 
potential 

-5e~ 5r 

V (r) = j—^, (21) 

when the screening parameter 6 > is not too large. To employ the FM we expand the 
Hulthén potential into the Laurent series 

V(r) = -^ + V reg (r), (22) 

with the regular part given by 

r oo 

Vr eg {r) = --5Y J 9nM 2n+ \ (23) 

Z ;i=0 



Application of the Fróbenius method 



8 



where 

^ [2(n + l)]!' 1 J 

and /3 n are given by a convenient expression 

2n-l i A: 

^=(-1)"^^ E •-')'(■':)/" '• (25) 




The convergence radius of (J2l?j) is py = -y, and in the range < r < py the potential 
(j22|) can be approximated with an arbitrary accuracy by a series with a finite number of 
terms. With the series truncated after the r 2P+1 term, the recursion relation (|12j) takes 
a form 

-2a,-! + 2(f - A)a t _ 2 - 2 ( . 

üi ~ i(i + 2l + l) • 

In Table |3]the numerical results are presented for different vàlues of 6 at various radii of 
enclosure R < py- Numerical stability was achieved by increasing the number of terms 
both in the series solution (JTTjl and in the potential expansion ()23|) until the approximate 
vàlues X n i for fixed R become stable to the quoted accuracy. The table also contains 
the exact eigenvalues E n i for the unconfined Hulthén potential, obtained analytically 
(l = 0) and by numerical integration (/ ^ 0) |15j . 



Table 3. Bound state energies of the Hulthén potential for various vàlues of the 
radius of confincmcnt R. 



1 


n = n r 


s 


A ni (4) 


A„i(6) 


A ni (8) 


A nZ (12) 


E n l 








0.050 


-0.4585448 


-0.4752873 


-0.4753117 


-0.4753125 


-0.4753125 






0.075 


-0.4463941 


-0.4631776 


-0.4632024 


-0.4632031 


-0.4632031 





1 


0.050 


0.4447886 


-0.0606327 


-0.0888523 


-0.0975630 


-0.1012500 






0.075 


0.4567303 


-0.0492481 


-0.0776572 


-0.0864942 


-0.0903125 


1 





0.050 


0.1680730 


-0.0802458 


-0.0947675 


-0.0992341 


-0.1010425 






0.075 


0.1800057 


-0.0687394 


-0.0834014 


-0.0879576 


-0.0898478 



4. Unconfined potentials 

Now we come to the discussion of unconfined potentials, namely to the case of a particle 
with the angular momentum / in a spherically symmetric potential of the form (@J) 
without any external enclosure. An unconfined system can be also effectively treated 
by the FM, as demonstrated by the calculation of the ground state energy of the one- 
dimensional anharmonic oscillator to an enormous precision of 1184-digits jH]. Here we 
show that the method can be implemented in a way that also allows for computing the 
excited states energies. We take the generalized series (fTT|). as a solution of the radial 
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Schròdinger equation, which fulfils the boundary condition u(X, 0) = 0, and consider two 
ways of imposing the second boundary condition at finite r = R, namely u(R, A) = or 
u^(R, À) = 0. At R — > oo both conditions are satisfied at the same vàlues of À, which 
correspond to the energy eigenvalues E n \ of the unconfined system, where n = 0,1,2, .... 
At finite R the vàlues of A are different: let us denote by A/w the vàlues satisfying 

u(R,X kl ) = for fe = 0,1,..., (27) 

and by X kl those satisfying 

uU(R,X' kl ) = 0, for fe = 0,1,... (28) 

In both cases we deal with the Sturm-Liouville eigenvalue problem; the nodes of the 
function u(r, Xm) in the radial variable divide thus the domain (0, R) precisely into fe- 
parts jUIj, and the same is true for the function u(r, X' kl ). The energy eigenvalues satisfy 
thus the following inequalities: 

X m < X i < X' u < Xn < X' 2l < X 2i < ... < X H < X k i < ••• (29) 

It is easy to check that the sign of u(R, X) is opposite to that of u^(R, X), for any value 
A lying within (X' kl ,Xki), which will be called the bounding interval in the following. 
Whereas, for vàlues of A in the interval (Xki, X k+li ) the signs of u(R, X) and u^'(R, X) 
are the same. These properties are best illustrated on FigJSl in the example of u(r, X) 
calculated for the spherical harmònic oscillator with angular momentum 1 = 0. 

1 

0.75 
0.5 
0.25 


-0.25 
-0.5 
-0.75 

2 4 6 8 10 12 

1 

Figure 3. Behavior of u(R,\) (sòlid line) and u^(R, X) (dashed line) for the 
spherically symmetric harmònic oscillator with u> = 1 and í = 0atü=2.5. The 
three first bounding intervals are marked by thick line 

Bound states in the radial potential V(r) are possible only for A < lim r ^ oo y(r). 
In this case, there exists a point R c such that for r > R c we have V e ff(r, l) — X > 
and the sign of u^ 2 \r, X) is the same as that of u(r, X), which implies that for r > R c 
the function u(r, X) must neither have a local maximum if u(R c , X) > 0, nor a local 
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minimum if u(R c , X) < 0. If both u(R c , A) and u^ l \R c , X) are positive (negative), then 
u(r, A) and u^'(r, A) tend monotonically to the plus (minus) infinity. Therefore, if the 
point of imposing the boundary condition, R is larger than R c , the following inequality 
is satisfied 

Ki < Eoi < ^01 < ■■■ < X kl < E kí < X k i < ... < \ ni < E n i < A„i, (30) 

where A„,/ < V e ff(R, ï). For increasing R, the energy X' kl (R) grows and Xki(R) decreases, 
approaching the exact eigenvalue from both sides monotonically. This allows for 
bounding the energy eigenvalues of an unconfined system with a required precision. With 
the approximate u(R, X), obtained by truncating the number of terms in the generalized 
power series to K, both pTj) and (}28|) are polynomial equations in A. The large set of 
energy levels can be thus determined numerically by finding the roots of polynomials. 
Generally, the bounding intervals are larger for higher states but for increasing R all 
the bounding interval shrink. We determine thus the new bounding intervals with the 
value of R increased by a suitably chosen Ai?. The iteration procedure is repeated for 
Ri = R + iAR, until the bounding energies for a chosen state (k, l) become equal to 
the accuracy desired, which determines its energy with that accuracy. If the procedure 
does not converge for the state of interest, the number of terms K, which are included 
in the power series (fTTj) should be increased. 



4-1. The Kratzer potentíal 

For demonstrating the convergence of the algorithm previously formulated, we first 
consider the Schròdinger problem in the Kratzer potential 

V (r) = + *=i, (31) 
for which the exact energy levels are given by 



E nl = -2d\(2n + 1 + ^{21 + l) 2 + 8cL 2 )- 2 , (n = 0, 1, 2, ...) (32) 

Choosing cLx = —8, and cL 2 = 4 as the parameters of the potential (jïïTJ), we carry out 
the calculation for the state (1,1), taking the number of terms in the power series (|TT|) 
suitably large (K = 160) in order to assure the numerical stability. The Table 0] shows 
with 14— digit precision the vàlues of Xn(Ri), obtained from the condition (|27j) . and 
those of A 11 (i?j), obtained from ([28)1 . for Rj = R + iAR. The bounding interval shrinks 
very fast and the exact value of energy is easily determined with 14-digit precision, 
En — —1.4476568219254. Figure 0] shows how the deviations from the exact energy 
monotonically diminish for increasing R, the energy difference Xu(R) — En approaches 
zero from above, and X'u(R) — En from below. 
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Table 4. Comparison of the convergence of Xn(R) and X 11 (R) to the exact value 



En = -1.4476568219254 for different vàlues of R and AR 



i R AR Au(i2i) A' n (Ei) 

5 0J5 -1.3178526137388 -1.8092879070838 

1 -1.3839260262988 -1.6834926644347 

2 -L4178923369387 -1.5810678698537 

3 -L4345 1425 11694 -1.5094745737014 

8 ï -1.4468612319721 -1.4500438665400 

1 - 1.447 5622942994 - 1.447 8789515475 

2 - 1.4476 472391241 - 1.4476 761151204 

3 - 1.44765 59628487 - 1.44765 83796562 

14 2 - 1.44765682 15565 - 1.44765682 24842" 

1 - 1.44765682192 38 - 1.44765682192 76 

2 -1.4476568219254 -1.4476568219254 




11.4 11.6 11.8 12 12.2 12.4 12.6 12.8 



Figure 4. Convergence of Àn(iZ) — 25ii, and X 11 (R) — En to zero, demonstrated in 
function of R 

4-2. Anharmonic oscillator 

The unconfined quantum anharmonic oscillator is one of the most frequently discussed 
quantum systems. Especially the one-dimensional example, being the simplest tractable 
model of quantum field theory [16J, is routinely used for examining the validity of various 
approximation methods ^7j, as its exact solution can be numerically determined to 
an arbitrary accuracy jHl HH1 EU 1201- The D— dimensional case is much less studied, 
generally, a spherically symmetric anharmonic potential ()15j) for different vàlues of the 



Application of the Fróbenius method 



12 



anharmonicity power J is discussed. The quartic potential (J = 2) was studied by 
means of various approximation methods, e.g. the self-similar approximation [21], the 
random phase approximation [22], and the artificial perturbation method [23] • There 
are also a few reports in the literature of the numerically exact results for the quartic 
[IBI 121 123, sextic (J = 3) [HEU, and octic (J = 4) potential [23 in the limited range 
of the anharmonicity parameter z. 

Here we show that the FM enables us to determine effectively the spectrum of the 
unconfined spherical anharmonic oscillator ()15j) in the wide range of the parameter z. 
With the solution of the Schròdinger equation u(r, À) in the form of the generalized 
series (119)1 we determine the numerical vàlues of bound-state energies, using the 
procedure formulated in the beginning of this section. After checking that the results 
available in the literature [H EU EH] are easily recovered to the quoted accuracy, we 
performed an extensive calculation of the spectrum of spherically symmetric anharmonic 
oscillators. The highly accurate results presented here may serve for testing the quality 
of various approximation methods. Especially challenging test is provided by the data 
obtained for negative vàlues of the parameter z, when the anharmonic potential has 
a Mexican hat shape. This range was not explored before, and the numerical data 
for bound-state energies were lacking. Therefore, in table 03 we present our results for 
several lowest states (n, /) energies of the quartic oscillator ( J = 2) at negative vàlues 
of z. In table E3 we compare the bound-state energies for oscillators with different 
anharmonicity powers (J = 2,3,4) at various vàlues of z. The results up to 8-decimal 
precision are quoted, but it is not difficult to improve the accuracy at will. However, one 
has to note that the efíbrt of the calculation increases for higher states. The appropriate 
shrinking of the intervals bounding the higher states energies is achieved only at large 
R, which requires the larger number of terms K to be included in the power series. The 
computational effort increases strongly, when the parameter z becomes more negative, 
i.e. for increasing radius of the hat. For example, equality of the vàlues Àoo(-R) and 
\' 00 (R) with 8-digit accuracy was achieved for z = —2 at the radius R = 3, which 
requires K = 65, while for z = —10 the same is obtained only at R = 3.9, which 
requires K = 140 . 

Table 5. The lowest bound state energies of the quartic oscillator for negative vàlues 



of z. 


Z 


-E"oo 










-10 


-21.88965823 


-21.66759959 


-21.22775192 


-20.57776933 


-19.72742347 


-9 


-17.30811508 


-17.05634313 


-16.56006280 


-15.83155272 


-14.88543607 


-8 


-13.23757818 


-12.94630409 


-12.37723745 


-11.55086571 


-10.48917435 


-7 


-9.68064055 


-9.33403278 


-8.66866023 


-7.71939851 


-6.51858369 


-6 


-6.64062824 


-6.21133359 


-5.41635974 


-4.31340468 


-2.94717508 
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Table 6. The low-lying eigenvalues of anharmonic oscillators for certain powers of 
anharmonicity 2 J at different vàlues of the parameter z 



J 


z 








E\o 


.È/03 


2 


-5 


-4.11911133 


-3.55958100 


-2.59137576 


-0.21067602 


-1.30027463 




-4 


-2.10349462 


"1 O A 1 Tr"ooo 

-1.34175838 


-0.15162331 


1.57335020 


1.36098787 




-3 


-0.54212526 


0.50094741 


1.95580851 


3.17712298 


3.71671698 




-2 


0.66142890 


2.04064501 


3.78779723 


4.66539082 


5.81395433 




0.1 


2.46463653 


4.58321826 


6.96529241 


7.45891176 


9.56314673 




0.5 


2.73789227 


4.99143053 


7.49177505 


7.94240398 


10.19704822 




1 


3.05794573 


5.47591999 


8.12171427 


8.52673739 


10.95972825 




2 


3.63948205 


6.37163094 


9.29940583 


9.63362791 


12.39687279 




5 


5.06966367 


8.63688366 


12.33685580 


12.55011513 


16.15843962 




8 


6.21722563 


10.49477725 


14.87064591 


15.02454165 


19.33921639 


3 


-8 


-5.00982691 


-4.08581639 


-2.49927339 


1.27019537 


-0.37612149 




-5 


-1.39861613 


-0.06070241 


1.93123653 


4.23192390 


4.43629185 




-2 


l. 24144988 


3.20845649 


5.72356878 


7.05803262 


8.68689446 




-i 

-1 


1.94504486 


4.15238570 


6.86055253 


7.97745578 


9.98954703 




0.1 


2.63985266 


5.11853012 


8.04543812 


8.96977843 


11.36185573 




1 


3.15630057 


5.85836853 


8.96728095 


9.76393591 


12.44002444 


4 


-8 


-3.20712029 


-1.83434272 


0.36915524 


4.03170882 


3.24522274 




-5 


-0.57239228 


1.19930711 


3.74932588 


6.24632114 


6.93207343 




-2 


1.55575737 


3.85587069 


6.82325601 


8.47500053 


10.35857448 




-1 


2.16525909 


4.66276445 


7.78250948 


9.21918447 


11.44445127 




0.1 


2.78594828 


5.50821786 


8.80159788 


10.03495445 


12.60737664 




1 


3.25887003 


6.16887293 


9.60819432 


10.69825337 


13.53482309 


5 


0.1 


2.91188167 


5.81747881 


9.38031919 


10.83477815 


13.54422429 




1 


3.35601445 


6.42909865 


10.11830492 


11.42342067 


14.38432780 



Our results for the quartic oscillator are presented graphically in FigEl The 
excitation energies with respect to the ground state, E n i — Eo Q , are plotted in function of 
the parameter z, which covers the range —15 < z < 25. The vàlues of z < correspond 
to the Mexican hat, and those of z > to the single-well shape of the anharmonic 
potential (|T7)l. It is interesting to note that, in spite of this difference, the dependence 
of excitation energies on z is smooth at the point separating the two cases, z — 0, which 
corresponds to the strong coupling limit, g — > oo. Instead, the behavior of excitation 
energies in the weak coupling limit (|^| — > oo) is very different for positive and negative 
z. For z— >oc we recover the radial harmònic oscillator of the frequency 2z, which energy 
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Figure 5. The behavior E n i — Eqq in function of z, the asymptotic behavior of 
E n — Eq — n^J—Az is shown by dashed lines. Two points of level crossing are marked 
by circles. 



eigenvalues are given by 

E nl — -> (2n + l + ^)V2~z. (33) 
z — > oo 

In FigEl one can observe how the levels with the same 2n + l become degenerate for 
increasing z. For z < the Mexican hat shaped anharmonic potential has a minimum 
at r min = J=f. If >> l 2 , the minimum of V e ff is close to r min and the effective 
potential is approximated well by 

K „ H ,4 + ^ )(r _^ f , (34) 

which does not depend on /, and corresponds to the one-dimensional harmònic oscillator 
of the frequency y/—4z. Therefore, in the limit 2^—00 we have 

E n — > + (n+^) v /z 4i(n = 0,l,2,...) (35) 

z — > —00 

In Fig|3]tlie grouping of states with different quantum numbers / but the same value of 
n can be observed for strongly negative vàlues of z. For z — > —00 all the states in the 
group approach the asymptotic behavior of E n — E = ri\f—Az. 

One can also notice an interesting phenomenon of level crossing that appears in 
the range of z < 0: at certain negative value of z two adjacent eigenvalues become 
degenerate. The two examples in Figj^l are marked by circles: the crossing point 
for the states (1,0), (0,3), which appears at z\ = —3.73656382, and the crossing 
point for (1,0), (0,4), which appears at z 2 = —5.42007803. The configuration of five 
lowest states, which for z > z\ is given by (0, 0), (0, 1), (0, 2), (1, 0), (0, 3), changes into 
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(0, 0), (0, 1), (0, 2), (0, 3), (1, 0) for z 2 < z < z x , and into (0, 0), (0, 1), (0, 2), (0, 3), (0, 4) 
for z < Z2- It should be stressed that in the case of one-dimensional anharmonic 
oscillator the phenomenon of level crossing does not appear, this becomes possible only 
for systems of D > 2 dimensions. 

5. Conclusion 

We have shown that the application of the Fróbenius method to the spherically 
symmetrical potentials of the form V(r) = + — + J2iZo allows an easy 
determination of the energy spectrum. This was demonstrated first for systems enclosed 
in a spherical box of the radius R, by studying the confined harmònic and anharmonic 
oscillators and Hulthén potential. With the increasing radius of confmement the bound 
states energies have been shown to approach those of the corresponding unconfmed 
systems. Even in the case of the Hulthén potential, when the convergence radius of the 
potential is finite, we obtain quite good approximations to the low-lying spectrum of the 
unconfmed potential, if the screening is not too strong. Later, we developed an efficient 
scheme for computing eigenvalues of the Schròdinger equation for unconfmed potentials 
with a controlled accuracy. The method allowed us to determine the low-lying states 
energies of spherically symmetric anharmonic oscillators with very high accuracy and 
moderate effort. Determination of energies becomes computationally more demanding 
for higher states, since more terms have to be included in the generalized power series. 
Our calculations cover a broad range of anharmonic parameters, both in the case of 
single well potential and for the Mexican hat shape. In the later case the computational 
effort increases strongly with the increasing radius of the hat. 

The method presented in this work can be easily applied for computing the precise 
spectrum of other spherically symmetric potentials. In the present work, we have studied 
the case of three-dimensional space but the calculation of energy levels in the two- 
dimensional case can be performed along similar lines. The results in even and odd 
higher dimension D can be easily derived from those in 2— and 3— dimensional space, 
respectively, via the transformation l — > 1+^=2. . One has to add that the method can be 
also used to determine the approximate wave functions. After substituting the calculated 
bound-state energy to the recursion relation (|T2~j) the coefficients of the generalized 
series (fTT|) can be successively determined in order to obtain the unnormalized wave 
function as a sum of the series. 
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